The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Geographic source of bats killed at wind-energy facilities in the eastern United States.

Bats subject to high rates of fatalities at wind-energy facilities are of high conservation concern due to the long-term, cumulative effects they have, but the impact on broader bat populations can be difficult to assess. One reason is the poor understanding of the geographic source of individual fatalities and whether they constitute migrants or more local individuals. Here, we used stable hydrogen isotopes, trace elements and species distribution models to determine the most likely summer geographic origins of three different bat species (Lasiurus borealis, L. cinereus, and Lasionycteris noctivagans) killed at wind-energy facilities in Ohio and Maryland in the eastern United States. In Ohio, 41.6%, 21.3%, 2.2% of all individuals of L. borealis, L. cinereus, and L. noctivagans, respectively, had evidence of movement. In contrast, in Maryland 77.3%, 37.1%, and 27.3% of these same species were classified as migrants. Our results suggest bats killed at a given wind facility are likely derived from migratory as well as resident populations. Finally, there is variation in the proportion of migrants killed between seasons for some species and evidence of philopatry to summer roosts. Overall, these results indicate that the impact of wind-energy facilities on bat populations occurs across a large geographic extent, with the proportion of migrants impacted likely to vary across species and sites. Similar studies should be conducted across a broader geographic scale to understand the impacts on bat populations from wind-energy facilities.

Genetic mechanisms and biological processes underlying host response to ophidiomycosis (snake fungal disease) inferred from tissue-specific transcriptome analyses.

Emerging infectious diseases in wildlife species caused by pathogenic fungi are of growing concern, yet crucial knowledge gaps remain for diseases with potentially large impacts. For example, there is detailed knowledge about host pathology and mechanisms underlying response for chytridiomycosis in amphibians and white-nose syndrome in bats, but such information is lacking for other more recently described fungal infections. One such disease is ophidiomycosis, caused by the fungus Ophidiomyces ophidiicola, which has been identified in many species of snakes, yet the biological mechanisms and molecular changes occurring during infection are unknown. To gain this information, we performed a controlled experimental infection in captive Prairie rattlesnakes (Crotalus viridis) with O. ophidiicola at two different temperatures: 20 and 26°C. We then compared liver, kidney, and skin transcriptomes to assess tissue-specific genetic responses to O. ophidiicola infection. Given previous histopathological studies and the fact that snakes are ectotherms, we expected highest fungal activity on skin and a significant impact of temperature on host response. Although we found fungal activity to be localized on skin, most of the differential gene expression occurred in internal tissues. Infected snakes at the lower temperature had the highest host mortality whereas two-thirds of the infected snakes at the higher temperature survived. Our results suggest that ophidiomycosis is likely a systemic disease with long-term effects on host response. Our analysis also identified candidate protein coding genes that are potentially involved in host response, providing genetic tools for studies of host response to ophidiomycosis in natural populations.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Low-coverage whole genome sequencing for highly accurate population assignment: Mapping migratory connectivity in the American Redstart (Setophaga ruticilla).

Understanding the geographic linkages among populations across the annual cycle is an essential component for understanding the ecology and evolution of migratory species and for facilitating their effective conservation. While genetic markers have been widely applied to describe migratory connections, the rapid development of new sequencing methods, such as low-coverage whole genome sequencing (lcWGS), provides new opportunities for improved estimates of migratory connectivity. Here, we use lcWGS to identify fine-scale population structure in a widespread songbird, the American Redstart (Setophaga ruticilla), and accurately assign individuals to genetically distinct breeding populations. Assignment of individuals from the nonbreeding range reveals population-specific patterns of varying migratory connectivity. By combining migratory connectivity results with demographic analysis of population abundance and trends, we consider full annual cycle conservation strategies for preserving numbers of individuals and genetic diversity. Notably, we highlight the importance of the Northern Temperate-Greater Antilles migratory population as containing the largest proportion of individuals in the species. Finally, we highlight valuable considerations for other population assignment studies aimed at using lcWGS. Our results have broad implications for improving our understanding of the ecology and evolution of migratory species through conservation genomics approaches.

Genetic characterization of potential venom resistance proteins in California ground squirrels (Otospermophilus beecheyi) using transcriptome analyses.

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

Inferring population connectivity in eastern massasauga rattlesnakes (Sistrurus catenatus) using landscape genetics.

Assessing the environmental factors that influence the ability of a threatened species to move through a landscape can be used to identify conservation actions that connect isolated populations. However, direct observations of species' movement are often limited, making the development of alternate approaches necessary. Here we use landscape genetic analyses to assess the impact of landscape features on the movement of individuals between local populations of a threatened snake, the eastern massasauga rattlesnake (Sistrurus catenatus). We linked connectivity data with habitat information from two landscapes of similar size: a large region of unfragmented habitat and a previously studied fragmented landscape consisting of isolated patches of habitat. We used this analysis to identify features of the landscape where modification or acquisition would enhance population connectivity in the fragmented region. We found evidence that current connectivity was impacted by both contemporary land-cover features, especially roads, and inherent landscape features such as elevation. Next, we derived estimates of expected movement ability using a recently developed pedigree-based approach and least-cost paths through the unfragmented landscape. We then used our pedigree and resistance map to estimate resistance polygons of the potential extent for S. catenatus movement in the fragmented landscape. These polygons identify possible sites for future corridors connecting currently isolated populations in this landscape by linking the impact of future habitat modification or land acquisition to dispersal ability in this species. Overall, our study shows how modeling landscape resistance across differently fragmented landscapes can identify habitat features that affect contemporary movement in threatened species in fragmented landscapes and how this information can be used to guide mitigation actions whose goal is to connect isolated populations.

Venom Gene Sequence Diversity and Expression Jointly Shape Diet Adaptation in Pitvipers.

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA2s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms-such as higher absolute levels of expression-are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Individual Variability in Bothropsatrox Snakes Collected from Different Habitats in the Brazilian Amazon: New Findings on Venom Composition and Functionality.

Differences in snake venom composition occur across all taxonomic levels and it has been argued that this variation represents an adaptation that has evolved to facilitate the capture and digestion of prey and evasion of predators. Bothrops atrox is a terrestrial pitviper that is distributed across the Amazon region, where it occupies different habitats. Using statistical analyses and functional assays that incorporate individual variation, we analyzed the individual venom variability in B. atrox snakes from four different habitats (forest, pasture, degraded area, and floodplain) in and around the Amazon River in Brazil. We observed venom differentiation between spatially distinct B. atrox individuals from the different habitats, with venom variation due to both common (high abundance) and rare (low abundance) proteins. Moreover, differences in the composition of the venoms resulted in individual variability in functionality and heterogeneity in the lethality to mammals and birds, particularly among the floodplain snakes. Taken together, the data obtained from individual venoms of B. atrox snakes, captured in different habitats from the Brazilian Amazon, support the hypothesis that the differential distribution of protein isoforms results in functional distinctiveness and the ability of snakes with different venoms to have variable toxic effects on different prey.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake.

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Predicting migration routes for three species of migratory bats using species distribution models.

Understanding seasonal variation in the distribution and movement patterns of migratory species is essential to monitoring and conservation efforts. While there are many species of migratory bats in North America, little is known about their seasonal movements. In terms of conservation, this is important because the bat fatalities from wind energy turbines are significant and may fluctuate seasonally. Here we describe seasonally resolved distributions for the three species that are most impacted by wind farms (Lasiurus borealis (eastern red bat), L. cinereus (hoary bat) and Lasionycteris noctivagans (silver-haired bat)) and use these distributions to infer their most likely migratory pathways. To accomplish this, we collected 2,880 occurrence points from the Global Biodiversity Information Facility over five decades in North America to model species distributions on a seasonal basis and used an ensemble approach for modeling distributions. This dataset included 1,129 data points for L. borealis, 917 for L. cinereus and 834 for L. noctivagans. The results suggest that all three species exhibit variation in distributions from north to south depending on season, with each species showing potential migratory pathways during the fall migration that follow linear features. Finally, we describe proposed migratory pathways for these three species that can be used to identify stop-over sites, assess small-scale migration and highlight areas that should be prioritized for actions to reduce the effects of wind farm mortality.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors.

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Pedigree-based assessment of recent population connectivity in a threatened rattlesnake.

Managing endangered species in fragmented landscapes requires estimating dispersal rates between populations over contemporary timescales. Here, we developed a new method for quantifying recent dispersal using genetic pedigree data for close and distant kin. Specifically, we describe an approach that infers missing shared ancestors between pairs of kin in habitat patches across a fragmented landscape. We then applied a stepping-stone model to assign unsampled individuals in the pedigree to probable locations based on minimizing the number of movements required to produce the observed locations in sampled kin pairs. Finally, we used all pairs of reconstructed parent-offspring sets to estimate dispersal rates between habitat patches under a Bayesian model. Our approach measures connectivity over the timescale represented by the small number of generations contained within the pedigree and so is appropriate for estimating the impacts of recent habitat changes due to human activity. We used our method to estimate recent movement between newly discovered populations of threatened Eastern Massasauga rattlesnakes (Sistrurus catenatus) using data from 2996 RAD-based genetic loci. Our pedigree analyses found no evidence for contemporary connectivity between five genetic groups, but, as validation of our approach, showed high dispersal rates between sample sites within a single genetic cluster. We conclude that these five genetic clusters of Eastern Massasauga rattlesnakes have small numbers of resident snakes and are demographically isolated conservation units. More broadly, our methodology can be widely applied to determine contemporary connectivity rates, independent of bias from shared genetic similarity due to ancestry that impacts other approaches.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complex-ity can be characterized by testing hypothesized links betweencomplex forms and their functions across species. Predatory ven-oms are composed of multiple proteins that collectively function toincapacitate prey. Venom complexity fluctuates over evolutionarytimescales, with apparent increases and decreases in complexity,and yet the causes of this variation are unclear. We tested alterna-tive hypotheses linking venom complexity and ecological sourcesof selection from diet in the largest clade of front-fanged ven-omous snakes in North America: the rattlesnakes, copperheads,cantils, and cottonmouths. We generated independent transcrip-tomic and proteomic measures of venom complexity and collatedseveral natural history studies to quantify dietary variation. Wethen constructed genome-scale phylogenies for these snakes forcomparative analyses. Strikingly, prey phylogenetic diversity wasmore strongly correlated to venom complexity than was overallprey species diversity, specifically implicating prey species’ diver-gence, rather than the number of lineages alone, in the evolutionof complexity. Prey phylogenetic diversity further predicted tran-scriptomic complexity of three of the four largest gene familiesin viper venom, showing that complexity evolution is a concertedresponse among many independent gene families. We suggest thatthe phylogenetic diversity of prey measures functionally relevantdivergence in the targets of venom, a claim supported by sequencediversity in the coagulation cascade targets of venom. Our resultssupport the general concept that the diversity of species in an eco-logical community is more important than their overall number indetermining evolutionary patterns in predator trait complexity.

Individual variability in bothrops atrox snakes collected from different habitats in the Brazilian Amazon: new findings on venom composition and functionality

Differences in snake venom composition occur across all taxonomic levels and it has been argued that this variation represents an adaptation that has evolved to facilitate the capture and digestion of prey and evasion of predators. Bothrops atrox is a terrestrial pitviper that is distributed across the Amazon region, where it occupies different habitats. Using statistical analyses and functional assays that incorporate individual variation, we analyzed the individual venom variability in B. atrox snakes from four different habitats (forest, pasture, degraded area, and floodplain) in and around the Amazon River in Brazil. We observed venom differentiation between spatially distinct B. atrox individuals from the different habitats, with venom variation due to both common (high abundance) and rare (low abundance) proteins. Moreover, differences in the composition of the venoms resulted in individual variability in functionality and heterogeneity in the lethality to mammals and birds, particularly among the floodplain snakes. Taken together, the data obtained from individual venoms of B. atrox snakes, captured in different habitats from the Brazilian Amazon, support the hypothesis that the differential distribution of protein isoforms results in functional distinctiveness and the ability of snakes with different venoms to have variable toxic effects on different prey.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option

Using trace elements to identify the geographic origin of migratory bats.

The expansion of the wind energy industry has had benefits in terms of increased renewable energy production but has also led to increased mortality of migratory bats due to interactions with wind turbines. A key question that could guide bat-related management activities is identifying the geographic origin of bats killed at wind-energy facilities. Generating this information requires developing new methods for identifying the geographic sources of individual bats. Here we explore the viability of assigning geographic origin using trace element analyses of fur to infer the summer molting location of eastern red bats (Lasiurus borealis). Our approach is based on the idea that the concentration of trace elements in bat fur is related through the food chain to the amount of trace elements present in the soil, which varies across large geographic scales. Specifically, we used inductively coupled plasma-mass spectrometry to determine the concentration of fourteen trace elements in fur of 126 known-origin eastern red bats to generate a basemap for assignment throughout the range of this species in eastern North America. We then compared this map to publicly available soil trace element concentrations for the U.S. and Canada, used a probabilistic framework to generate likelihood-of-origin maps for each bat, and assessed how well trace element profiles predicted the origins of these individuals. Overall, our results suggest that trace elements allow successful assignment of individual bats 80% of the time while reducing probable locations in half. Our study supports the use of trace elements to identify the geographic origin of eastern red and perhaps other migratory bats, particularly when combined with data from other biomarkers such as genetic and stable isotope data.

Venom variation in Bothrops asper lineages from North-Western South America.

Bothrops asper is a venomous pitviper that is widely distributed and of clinical importance in Mesoamerica and northern South America, where it is responsible for 50-80% of all envenomations by Viperidae species. Previous work suggests that B. asper has a complex phylogeographic structure, with the existence of multiple evolutionarily distinct lineages, particularly in the inter-Andean valleys of north South America. To explore the impact of the evolutionary history of B. asper on venom composition, we have investigated geographic variation in the venom proteome of this species from the populations from the Pacific side of Ecuador and south-western Colombia. Among the 21 classes of venom components identified, proteins from mainly four major toxin families, snake venom metalloproteases (PI- and PII-SVMP), phospholipases A2 (K49- and D49-PLA2s), serine proteinases (SVSP), and C-type lectins-like (CTL) proteins are major contributors to the geographic variability in venom. Principal component analyses demonstrate significant differences in venom composition between B. asper lineages previously identified through combination of molecular, morphological and geographical data, and provide additional insights into the selection pressures modulating venom phenotypes on a geographic scale. In particular, altitudinal zonation within the Andean mountain range stands out as a key ecological factor promoting diversification in venom. In addition, the pattern of distribution of PLA2 molecules among B. asper venoms complements phylogenetic analysis in the reconstruction of the dispersal events that account for the current biogeographic distribution of the present-day species' phylogroups. Ontogenic variation was also evident among venoms from some Ecuadorian lineages, although this age-related variation was less extreme than reported in B. asper venoms from Costa Rica. The results of our study demonstrate a significant impact of phylogenetic history on venom composition in a pitviper and show how analyses of this variation can illuminate the timing of the cladogenesis and ecological events that shaped the current distribution of B. asper lineages. BIOLOGICAL SIGNIFICANCE: Bothrops asper, called "the ultimate pitviper" due to its defensive behavior, large body size, and medical importance, represents a species complex that is widely distributed from southern México southwards across north-western South America to north-western Perú. This work reports the characterization of the venom proteomes of B. asper lineages from the Pacific sides of Ecuador and south-western Colombia. Multivariate analyses indicate that variability in venom composition among the B. asper lineages is driven by proteins from four major toxin families, presumably in response to selection pressures created by recent and historical ecological conditions created by geological and climatic events from the Pliocene-Pleistocene to the present along the Central and South American Continental Divide. The emerging biogeographic pattern of venom variation, interpreted in the context of the current phylogenetic hypotheses, support and complement previously proposed evolutionary Plio-Pleistocene dispersal events that shaped the present-day distribution range of B. asper lineages. In addition, our venomics data indicate the occurrence of genetic exchange between Colombian and Pacific Costa Rican populations, which may have occurred during the second wave of B. asper migration into Mesoamerica. Our work represents a foundation for a future broader sampling and more complete "-omics" analyses to deepen our understanding of the patterns and causes of venom variation in this medically important pitviper.